Using Josephson junctions to determine the pairing state 
of superconductors without crystal inversion symmetry 
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Theoretical studies of a planar tunnel junction between two superconductors with antisymmet- 
ric spin-orbit coupling are presented. The half-space Green's function for such a superconductor 
is determined. This is then used to derive expressions for the dissipative current and the Joseph- 
son current of the junction. Numerical results are presented in the case of the Rashba spin-orbit 
coupling, relevant to the much studied compound CePtaSi. Current-voltage diagrams, differential 
conductance and the critical Josephson current are presented for different crystallographic orien- 
tations and different weights of singlet and triplet components of the pairing state. The main 
conclusion is that Josephson junctions with different crystallographic orientations may provide a 
direct connection between unconventional pairing in superconductors of this kind and the absence 
of inversion symmetry in the crystal. 

PACS numbers: 



I. INTRODUCTION 

The question of how parity violation affects supercon- 
ductivity has until recently not been subject to much ex- 
perimental studies. In recent years, however, supercon- 
ductivity has been discovered in several materials with 
a noncentrosymmetric crystal structure. This offers an 
arena for the study of superconductivity in the absence 
of inversion symmetry. Theoretical studies of such sys- 
tems have predicted several exotic features, reviewed in 
Rcfs. [HandS The absence of inversion symmetry allows 
an antisymmetric spin-orbit coupling in the Hamiltonian. 
This has, among other things, the consequence that the 
pairing state of the superconductor may not be classified 
as a spin singlet or a spin triplet state;" 

The most famous and studied example of the non- 
centrosymmetric superconductors is the heavy fermion 
compound CePtsSi, which possess several interest- 
ing properties ,ii2i5iS For instance, the pairing state of 
CePt3Si seems to contain line nodes^ even though NMR 
measurements are of the kind expected for a conventional 
superconductor.— Several theories have been put forward 
to explain thisi 7 -' 8 -' 9 - Other examples of noncentrosym- 
metric superconductors are Ulr, L^PdaB, I^PtaB, 
Cd2Re207 and possibly KOS2O6. The absence of inver- 
sion symmetry in these materials destroys spin degen- 
eracy through antisymmetric spin-orbit coupling. This 
is expected to be strong in some of the materials men- 
tioned above ) 10 ! 11 especially in compounds containing 
atoms with a large atomic number. Line nodes also 
seem to appear in the pairing state of I^PtsB ; 12 ' 13 
whereas Cd2Re2C>7, Li 2 Pd 3 B and KOs 2 06 appear to be 
nodeless . 12 ' 14 ' 15 

The experiments performed on these materials so far 
mostly concern quantities such as specific heat, magnetic 
penetration depth and the nuclear spin-lattice relaxation 
rate. They are all important in order to determine the 
pairing state of a superconductor. However, tunneling 
spectroscopy and experiments on Josephson junctions are 



also a very useful tool in this respect, both in conven- 
tional and high-T c superconductors . 16 ' 17 ' 18 

Recently, theoretical studies of transport in a junction 
between a normal metal and a noncentrosymmetric su- 
perconductor were performed ] 19 ' 20 These kind of trans- 
port measurements do not probe bulk properties directly, 
but will depend on how the pairing state is affected by the 
surface. Due to the possible triplet component, one may 
expect the gap to deviate from its bulk valu o 21 ' 22 and 
formation of Andreev bound states near the surface i 17 ' 23 

The transport properties of a Josephson junction con- 
sisting of two noncentrosymmetric superconductors have 
been investigated in Ref. [24]- Given particular pairing 
states, it was noted that both the quasiparticle current 
and the critical Josephson current would depend on the 
relative crystal orientation of the superconductors. Sim- 
ilar effects may appear with the two-band superconduc- 
tor MgB 2 j^' 2 £ However, in Ref. [H, the bulk density of 
states was used, neglecting the effect of surface scatter- 
ing. One might therefore question the validity of these 
results, since the effect of surface reflection was not con- 
sidered. 

In this paper, the effect of surface scattering is taken 
into account when determining the transport properties 
of the above mentioned Josephson junction. It is shown 
that the effects predicted in Ref. |2J may still appear, 
even though the surface provides a strong coupling be- 
tween the spin-orbit split bands. Thus, in some cases one 
may expect qualitative changes in the differential conduc- 
tance for different relative crystal orientations of the two 
superconductors. In addition, quantitative changes in the 
critical Josephson current may be expected. This could 
make it easier to establish a direct connection between 
the unconventional pairing and the absence of inversion 
symmetry. This paper is hence an attempt to motivate 
experimental work on such junctions. 

The paper is organized as follows. In section [Til we de- 
fine the model containing a general antisymmetric spin- 
orbit coupling. The Green's function is then established, 
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first in the bulk case and then in a half-space or semi- 
infinite scenario. Expressions for the tunneling currents 
are presented in section IIIII In section IIV| numerical re- 
sults using the Rashba spin-orbit coupling are presented 
to exemplify the predicted effects. 



II. 



MODEL 



We will start by considering the bulk properties of a 
clean superconductor with spin-orbit split bands. The 
model will be written down in the continuum limit. Hav- 
ing established the bulk Green's function, we move on to 
derive the Green's function in the presence of a reflecting 
surface. 



A. Bulk properties 

Let the Hamiltonian consist of two terms, H = Htsi + 
Hgc, a normal part and a part describing superconduc- 
tivity. In the bulk, the normal part is 



His 



dfc 



i [(en 



fj,)l + B k ■ a] <j> k , 



(1) 



The standard mean field approach gives 



H 



sc 



A 



dfc ( A\,kc{ k c\ k + A* X h cx,-hC\,h 



(3) 

where A x , k = -J2^Jdk'V\ IJ ,(k,k')(c^ k 'C. ti - k '). 

A\,—k = ~ A\ k follows from the fermionic anticom- 
mutation relations. One should note that the two 
bands are decoupled in the mean field approximation. 
However, the gaps A±, k are in general not independent, 
but related through the self-consistency equations due 
to the above mentioned possibility of interband pair 
scattering.— 

Let K. denote the time-reversal operator, whose effect 
on the operators in the spin basis is K : c k a — —&c^_ k _ a . 



It may be derived that K, : 



X.k 



*A,feC A ,-fe' where 



t\.k = —t\,-k is a gauge-dependent phase factor. One 
may write A Ajfe = iA,feXA,fe, where XA.fc is the order pa- 
rameter for pairs of time reversed states on which observ- 
able quantities will depend. Thus, x\,k = X\,~k may be 
expanded in terms of even basis functions of irreducible 
representations of the space group.— 

Define the matrix A k whose elements arc the gap func- 
tions Afc crcr' in a spin basis. By transforming equation 
([3]), one arrives at 



where </>£ = (eL., cL), s k is the band dispersion and /j, 
the chemical potential. The vector cr consists of the three 
Pauli matrices. 

The vector B k describes the anti-symmetric spin-orbit 
coupling. It removes the spin degeneracy from the band 
e k . The absence of inversion symmetry is reflected in 
the property B k = B k . An electron in a state with 
momentum fc will align its spin parallel or antiparallcl 
to B k . The symmetries of B k may be determined from 
point group symmetry considerations^ 

Diagonalization of equation ([1]) gives i?N = 
Y,\=±, k ^,kc\ k c x ,k, where £± ik = e k - ju ± \B k \. The 
spin of an electron in a state with momentum fc will point 
parallel(antiparallel) to B k in band +(— )■ 

We write down the term responsible for superconduc- 
tivity in terms of the long-lived excitations in the normal 
state, i.e. 

Hsc = \Y1 j dfcdfc'y AM (fc,fc') c A -fe4,fe^.fe'^-fe'- 

We will consider the limit where the spin-orbit splitting 
is much larger than the superconducting gaps. This is a 
relevant limit, at least for the materials CePtsSi^ and 
Cd2Re207»ii In that case, interband Cooper pairs are 
strongly suppressed, even though the two bands may 
touch at some isolated points on the Fermi surfaced 
Thus, the model Q contains only intraband Cooper pair- 
ing. However, it does include an internal Josephson cou- 
pling, i.e. scattering of Cooper pairs between the bands. 



A fc = ?7fc,s (-icy) + T] k ^ T (B k ■ cr) (-ia y ). 



(4) 



Thus, in the absence of spatial inversion symmetry, the 
order parameter in a spin basis have no definite parity, 
but is in general a linear combination of a singlet (S) and 
a triplet (T) part .JLi&Si The singlet and triplet compo- 
nents are determined by 



m,s = ijj (x+,h + x-.k) , 

Vk,T = ^ {X+,k - X-,k) ■ 



(5) 



There is no need to specify the momentum dependence 
of the gaps xx.k at this point. 

In the bulk, the Green's functions are diago- 
nal in momentum space due to translational sym- 
metry. In the imaginary time formalism, de- 
fine the normal and anomalous Green's functions as 
Gb^a'(k,r) = -(T T c kcr (T)c ka ,(0)) and T h ^ al {k,r) = 
(T T Cfeo-(r) c_fe CT '(0)), respectively, where the subscript b 
denotes bulk. It is convenient to transform to fermionic 
Matsubara frequencies, io n = (2n + l)7r//3, where (3 is the 
inverse temperature. The bulk Green's function in spin 
x particle-hole space, 

/ Q h (k,iu! n ) -T h (k,\uj n ) \ 
G b (fc,ko„) = . , (6) 

V iw n ) -gi(-k,-ioj n ) J 

is found by solving the Gor'kov equations (|A1[) presented 
in appendix The components are matrices in spin 
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space, given by 

Gb(k,iu) n ) = | ^2 c r ^ fc G A (A;,icj„) 



T h (k,iw n ) = -- ^ °% k a v F\(k,iuj n ) 



(7) 



A=± 



in terms of the complex scalar functions 



G\(k,iuj n ) = - 
F A (fc,iw n ) 
and the matrices 



XA.fc 

^n+Clfc + IXA.fcl 2 ' 



a% =t + \B k <r. 



(8) 



(9) 



B. Half-space Green's function 

The bands + and — defined in the previous section 
has the property that reversing the direction of an elec- 
tron's momentum while preserving its spin requires a 
change of bands. Thus, one would expect that the in- 
dependence of the bands + and — could be vulnerable 
to scattering, e.g. from impurities. In fact, it has been 
shown that a small concentration of nonmagnetic impu- 
rities does not change the picture of independent bands 
in the mean field approximation^ A perfectly reflect- 
ing surface should however lead to a severe mixing of 
the bands. This needs to be taken into account when 
describing transport in heterostructures containing these 
materials. 

The presence of a surface will make the Hamilto- 
nian and the Green's function nondiagonal in momentum 
space. Still, due to the nature of the spin-orbit coupling, 
it is convenient to work in a plane wave basis. We will 
assume that the surface is perfectly smooth. Of course, 
any real surface will have some roughness, which may 



J 



very well modify the results of this paper. However, at 
least for not too rough surfaces, this model is an ap- 
propriate starting point. We will also assume that the 
surface is spin inactive, i.e nonmagnetic. 

Consider the simplest case of a perfectly smooth sur- 
face at x = 0, such that the electrons are confined to 
x < 0. We seek the Green's function G(ki,k2,r), whose 
elements are 

g aa , (k u k 2 ,r) = -<T T c fel , CT (r)4 2j(T ,(0)) , 
^(fci.fa.r) = (T T c kl ,a{r)c-k 2 ,a'(0)) , (10) 

where Cfc liCr is the annihilation operator for a plane wave 
state. In the presence of a scattering surface, these corre- 
lation functions will not be diagonal in momentum space. 

Due to translational invariance in the y- and z- 
direction, it is natural to introduce the 4x4 Green's 
function in spin x particle-hole space in a mixed rep- 
resentation, G(xi,X2,ku,iu) n ). We have defined fey = 
k y y + k z z. The normal and anomalous components are 
<3**>(xi,X2,k h T) = -(T T c Xl> fe ||!(x (r)4 2ife||iCT ,(0)) and 

F<T<r'(xi,X2,k\\,T) = (T T c Xl ,k ll ,A T ) c z 2 ,-fcy ,<x' (°)) > re- 
spectively. The Green's function is determined by the 
Gor'kov equations (|A3|) . which are presented in appendix 
|A"1 The boundary conditions are 



, Xi — or x 2 = 0. 



(11) 



The pair potential in this mixed basis, A(xi, x 2 , fell), 
should be determined self-consistently. Even though it 
may deviate significantly from the bulk near surfaces ; 21 ! 22 
we will approximate it by its bulk value. This approxima- 
tion is expected to give qualitatively correct resultsj 17 i 20 
In appendix [X] it is shown that this approximation en- 
ables us to express the half-space Green's function in 
terms of bulk Green's functions. This may be realized 
by treating the surface as a wall of nonmagnetic impuri- 
ties of infinite strength.— The momentum space Green's 
function then becomes 



G(fci,fe 2 ,iw n ) = G b (fci,i^„)5(fci :X - /c 2 ,x) - G b (fci,i^„)G b 1 (0,0, fc||,iw„)Gb(A;2,iw„) S(k 1 ^ - k 2 ,\\). (12) 

I 



To determine this, we need the inverse of the matrix 

/DO 
dfc x G b (fc>„). (13) 
-oo 

Let us now define k = (k x ,k^) and k = (— k x ,ku). 
From the previous subsection, we saw that the gaps X±,k 
where unchanged upon reversal of the momentum. At 
this point we restrict ourselves to surfaces such that the 



gaps are unchanged also when reversing the component 
of the momentum perpendicular to the surface only, i.e. 
X\k = XA.fc- Although this is not a necessary require- 
ment to determine the Green's function, it will simplify 
the calculations and be sufficient for the scenarios consid- 
ered here. We will also assume £ A ^ = ^A.fe- Using these 
approximations, the properties G\(k,iuj n ) — G\(k,iu> n ) 
and F\(k,iw n ) — F\(k,iuj n ) follow from equations ©. 
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We now convert the fc^-integral in (|13[) to an energy in- 
tegral. The integrand will be strongly peaked about the 
Fermi level. Thus, we apply the quasiclassical approxi- 
mation of replacing all momentum-dependent quantities 
by their value at the Fermi levels We introduce the no- 
tation 



where N£ fep is |9fe x ^A,fc| 1 taken at k F . The vector 



Br 



(17) 



k F = (k F , x ,ku) , k F = (-k FtX ,k 



(14) 



where k F ^ x > is determined by £fc F = given fey. We 
define the quasiclassical or ^-integrated Green's functions 
by 



has the property b—j,.. = — , but is not a unit vec- 
tor. Similarly, the integral over the anomalous Green's 
function is obtained from 



gx(k F ,iu n ) 
f\(k F ,iw n ) 



\X\k F 



XA,fe F 



V^I + \x\k F \ 2 



(15) 



The integral over the normal Green's function in matrix 
(fl~3|) is found using 



1 f°° 

- J dk x a^ k G x (k, iuj n ) = irN^a^ g x (k F , \w r , 



(16) 



1 f°° 

- J dk x a% h F x (k,iu n ) = TrNl kw ag H fx(k F ,iuj n ). 

(18) 

We now assume that the difference in the density of states 
of the two spin-orbit split bands is small and may be 
neglected. Consequently, we also let N* kp = JVf kp = 
N£ f .— This is not a necessary step in order to proceed, 
but it simplifies the calculations. 



The inverse of the matrix (fT5)) is then 



a L ffp( fe F,iw„) \o p hk a y f p (k F ,icj n ) 



- ia y a b h „ /p( fc F,iw 



a v a b h v y g P (k F ,iuJn) 

re II 

^ + Re(x+,fc F X-, fcF ^ 



(19) 



(20) 



We have introduced the function 

K(k F ,iu n ) = 2 

where b± t h« = 1 db |&fe,, | 2 . Later, it will be apparent that zeros in K(k F ,iw n ) will correspond to surface bound states. 

Introduce the simplified notation G x ,i = G\(ki,uo n ), F p2 = F p ,(k2,uo n ) and g p = g p (k F ,iu) n ). No momentum 
index is needed on the latter since it depends only on the parallel momentum and fei ii = k 2i \\ = fc||. We are then 
ready to write down the half-space Green's function. The normal and anomalous parts, defined in (|10[) . are 



g(ki,k 2 ,iu n ) = ~ j ^CT^ fc GA,i<5(fci, K 



k 2 ,x) 



(21) 



\ V a Xf "* 

2nK(k F M n )N- F k ^ 



Gx 



Mg p .2 + f P F; ft ) + FxM;g»,2 - 9 P F; t2 )] W,, 



k 



and 
F{k 1 ,k 2 ,iuJ n ) 



-\^Y,^v F ^{K x -k 2 ^ x ) (2 



respectively. These functions are found by inserting equations ^ and (119|) into equation (|12D . 
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We have defined the matrix 



- 4 fcl <„ 4 fc2 - CI 1 + <m ■ - < ( 23 ) 



where the expressions for the scalar /3 k pfl f s)2 and the vector 
atjf^ are given in appendix [Bj 



III. CALCULATION OF TUNNELING 
CURRENTS 

Let us now consider a planar tunnel junction between 
two superconductors with spin-orbit split bands. We 
name the systems A and B and let the x-axis point per- 
pendicular to the junction. In addition, we use the letter 
c for operators and k for momenta on side A, and d and 
p for the corresponding quantities on side B. The spin- 
orbit coupling is described by the vectors B k and B p 
on each side. These vectors are not necessarily equal. 
Let us briefly exemplify this by considering the Rashba 
spin-orbit coupling, B k = a(h x fe), even though we will 
work with a general Bk. Here, the vector n describes 
the direction of broken inversion symmetry of the crys- 
tal. This means that if the crystallographic orientation 
on side B is different from side A, B^ and Bp will point 
in different directions even when k = p. 

The tunneling process is described by 



Cfcc 



Ht = Y.J dkd P { r k P ,aa'cl a d pal + T k * p aa ,dl a 

<J<j' 

. (24) 

The validity of results using perturbation theory in the 
tunneling Hamiltonian formalism has been showed by 
Pranged 

We emphasize that the systems are described in terms 
of plane wave states. Thus, Tk PtC!C! i is the transfer am- 
plitude from an incoming plane wave state with momen- 
tum p on side B to an outgoing plane wave state with 
momentum k on side A. When scattering a plane wave 
on a barrier, the perpendicular momentum of the trans- 
mitted wave points in the same direction as the incoming 
wave. In addition, we assume that the tunneling process 
conserves spin. These properties result in 



(25) 



where Q[x] is the Heaviside step function. Time- reversal 
symmetry also demands T* k p = Tfc iP . 

Of course, there is also an amplitude for the incom- 
ing plane wave being reflected. However, when treating 
equation (|24|) as a perturbation, the current will be ex- 
pressed as Green's functions of the unperturbed systems 
A and B. Thus, the reflection is taken into account by 
using the half-space Green's functions obtained in the 
previous section. 

The current from side B to side A is defined as J(i) = 
— e{AfA), where A/a is the total charge operator on side A 
and the operator A/a is given by the Heisenberg equation 



A/a = ^[Ht,Ma\- Treating the tunneling Hamiltonian 
as a perturbation, the Kubo formula gives I(t) — I qp + 
where 



J qp = -2eIm$(eV), 
7j(i)=2elm [e- 2ieVt V{eV)] 



(26) 



In the imaginary time formalism, when defining M(r) 
J2aa> / dfed P c L( T ) c W'0~): we have 



$(iw„) = - / dre iw " T (T r M(r)Mt(0)) , 
Jo 

f ■ 

*(iw„)=-/ dTe^ T (T T M(T)M(0)) . (27) 
Jo 



The time dependence of the operators are given by the 
unperturbed Hamiltonian and the expectation values are 
to be taken in the unperturbed state. The voltage 
is defined by eV = /ia — Mb- The bosonic Matsub- 
ara frequency is ui v = 2^7r//3, which will be subject to 
iu v -> eV + i0+. 

Equations (|27|) may be written as 



$(iw„) 



1 



dfeidfc 2 dpidp 2 T kuPl Tl 2 



(3 J "'" 1 "'" Z "- f ' ± " f ' i ~ re l>J'l~«2,JJ2 

x y^Tr^A(fc2,fci,i^n - iw„)0B(Pi,P2,iw n ) 



(28) 



and 



= iy dk 1 dk 2 dp 1 dp 2 T kuPl T^ p2 (29) 

,iw„ - \LO v )F B (p 1: p 2: \LLl n ) 



x E Tr 



where the components of the Green's functions are de- 
fined in ([TO]) and Tr denotes a trace in spin space. 



As before, when converting the momentum integrals 
to energy integrals, we replace all momentum-dependent 
quantities by their value on the Fermi level. The density 
of states N(kp) at the Fermi level is assumed equal in 
both bands. In addition, we assume N(kp) — N(kp) 



and T; 



k F ,pF 
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A. The quasiparticle current 



Inserting the Green's function (f2lj) in equation ([28)1 . 
one arrives at 



r 2 r' 



(30) 



*(iw v ) = ^ J dkj dp\T kFiPF \ 2 N A (k F )N B (p F ) 
i4^(fc F , PF )5 3 A ^(fc F ,p F ,i^) 

+ 1^^" (fe F , p F )Sl P ™ V (k F , PF , \W V ) 



where repeated greek indices are to be summed over. The 
prime indicates that the integrals over the Fermi surfaces 
are restricted to positive k x ,p x . The Ai 's are defined by^, 



A^(k,p) = Tra% A al B + Tra^JL , 

k p k V 

A^{k,p) = Tro^v™ + Tra% A a^ , 

k k 

A^(k,p) = Tr a^aZ B + Tra^a^ , (31) 

V V 

A x /^(k,p) = Tra^ffZ + ^5tT°lZ 



^ k.k P^P ^ k,k P<P 



These quantities depend on B^,B^,B?,Bp and ex- 
plicit expressions are given in appendix iBl The Si's de- 
pend on the momenta through the gaps and are defined 

as 

S^(k,p,uj v ) = i^^(fe,iw n -vu v )g*(p,iu> n ), 
S^ 17,u {k,p,iu) v ) = ~^2g£(k,iu n - iu v )T* vv (p, iw„) , 
S^ikypAuu) = -^r^(fc,iw„ - k^,).g^(p, iw n ) , 



-iXpfi'yrjv 



(k,p,iw v 



(32) 



The function iw„) was defined in equation (|15p 

The function T \ pPi {k,iui n ) is 



r ApM (fe,icj„) = 



5a [ffpSV + /p/;] + fx [fp9^ - 9pf£\ 



K(kAu n ) 



(33) 



where the arguments of the g's and f's were omitted for 
clarity. Note that r ApAt (fc F , iu> n ) does not depend on N% . 

In equation (|30|) . we have reached the point at which 
the current J qp is expressed as two surface integrals over 
half of the Fermi surface on each side. In addition, one 
is left with the Matsubara sums which may be converted 
to energy integrals. To get further, one must insert the 
appropriate angular dependence of the quantities x\,k F , 



-Bfep, N(k F ) on each side as well as |Xfc F 



In most 



cases, the remaining integrals need to be performed nu- 
merically. Both the energy and angle integrands contain 
integrable singularities which must be handled with care. 



We will now assume that the two gaps x±,k are phase- 
locked due to the internal Josephson coupling. We write 
out the phase explicitly, such that X±k ~~ * X± k e ^ ■ 
fc and X- k are rea ^ f rom now on > but not necessarily 
of the same sign. Obviously, the same also applies to the 
gaps on side B. 

To obtain the current / qp , we need lra^(iuj^). Since 
the Ai's are real (see appcndix lBj) . the only complex parts 
are contained in the Matsubara sums. By converting the 
sums to contour integrals in the complex plane and de- 
forming the contour, one finds that Im Si(eV + i0 + ) may 
be expressed as energy integrals containing the functions 
lmg x (k,E + i0+), ImT \ pfl (k, E + i0+) and the Fermi- 
Dirac distribution n F (E). Details of this procedure and 
the choice of appropriate branch cuts are found in ap- 
pendix [UJ The first function is proportional to the usual 
bulk density of states 



Im. 9A (fc,S + i0 4 



-e 



LEI 



\X\,k\ 



\E\ 



\X\,k\ 



(34) 

As before, Q[x] is the Heaviside step function. The sec- 
ond function is somewhat complicated, but may be writ- 
ten as 



i m r Aw (fe,E + io + ) = e 



LEI- 



Xm,fc | 



\Xm,k | 



E\ 6 - X+ ,kX-M Px m {E) S(\E\ - E Q t) , (35) 



where |x m ,fc| = min(|x +! fc|, X-,fc|)- The functions P\ PP (E) and P\ PP (E) are even functions of E and may be found 
by using equations (|C2|) and (IC3[) in appendix [Cl If one interprets Imr ApM (fe, E + i0 + ) as a density of states, the first 
term describes a continuum above the smallest gap. However, the second term describes additional discrete states 
below the smallest gap. These are the Andreev bound states induced by the reflection from the surface. Note that 
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they appear only when the sign of the two gaps differ, as was also noted in Ref. The energy E Q k is the positive 
solution to the equation 

b- M (x + MX-M - El k ) + 6 + , fc „ ^\X + A 2 -El^\x-A 2 -El k = (36) 



and is measured relative to the Fermi level. Thus, we 
get a band of low energy surface bound states in the part 
of momentum space where X+,kX—,k < 0- Equation (|36p 
corresponds to equation (14) of Ref. but here we have 
made no assumption of the particular form of . 

Let us also comment on what happens in the limit 
of a singlet superconductor. From equation l[5|). we 
see that this limit corresponds to x+.fe = X-M- ln 
that case, there are obviously no Andreev bound states 
and T\ ppi (k,\uj n ) = g + (k,\uj n )/A = g_(fe,iw„)/4. The 
current J qp then equals the result obtained using bulk 
Green's functioni 35 ' 36 

Whereas the limit x+,k = X-,k corresponds to a singlet 
superconductor, setting the gaps to zero corresponds to a 
normal metal. The current will in those cases not depend 
on the nature of the spin-orbit coupling vector The 
reader may wonder why there are no remnants of the 
spin-orbit coupling in these limits. This is a consequence 
of the approximation of equal densities of states for the 
two bands. 



B. The Josephson current 



The two-particle current is found by inserting the 
Green's function (f2"2")) in equation ([29]) . giving 



*(i^H— / dfc / dp\T kF , PF \ 2 N A (k F )N B (p F ) 



Al 1 (fe F , Pf ) <§i 7 ( &f , Pf , w u ) 



(37) 



~ ±A^(k F ,p F )S^(k F ,p F ,i^) 
+ i A x /^(k F ,p F )S^^(k F ,p F ,icj u ) 



As before, the integrals are restricted to positive k xl p x 
and repeated greek indices are summed over. The <SVs 



are defined as 

S^(k,p,iu u ) = i^/; A (fc,iw„ -bj u )f*(p,vu n ), 
S?f vv {k,p,iu v ) = i^/^ A (fc,iw„ - iui v )A B vv (p,iuj n ) . 

s^(k, P , = ~ A *^ k > iw « " ^ A B (P' . 



5' 



Xpfi-yrjiy 



(fc,p,iw„) 



(38) 







The function f\(k,\uj n ) is defined in (|15p and 
A ApA1 (fc,iw„) is 



K X p P {k,iui ri 



K(k,iui n ) 



(39) 

As in the previous subsection, we assume that the gaps 
are phase-locked, i.e. X±,k ~ * X±,k e ^ A an d treat X±h(v) 
as real. 

In the limit of a singlet superconductor, A\ pfJ _(k, iu n ) = 
/+(fc, iw„)/4 = /_(fe,iw n )/4. The Josephson current re- 
duces to the result found using bulk Green's functions* 3 -! 

It should be noted that (|37[) is a tunneling limit expres- 
sion. Thus, it may not capture all the unusual phenom- 
ena that arise when Andreev bound states contribute to 
Josephson currents^ 



IV. RESULTS 

In this section, we consider a junction consisting of two 
equal superconductors and present numerical results on 
the quasiparticle and Josephson current. We choose to 
study the Rashba interaction 



Bk = a(n x k) , 



(40) 



both because of its simplicity and its relevance to real ma- 
terials like CePt 3 Si (Ref. and Cd 2 Re 2 7 (Ref. [H). 
The vector h represents the direction of broken inversion 
symmetry of the crystal. 

We restrict ourselves to junctions where h A and h B are 
perpendicular to the tunneling direction, i.e. 
0. The angle £ is defined by 



cos C = n A • n B . 



(41) 
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Of course, from an experimental point of view, only dis- 
crete values of the angle £ may be realizable. 

The variation of the current with £ is a result of the 
facts that n determines the spin structure of the spin- 
orbit split bands and that spin is conserved in the tun- 
neling process. It should be noted that replacing one 
of the superconductors by a ferromagnet with magneti- 
zation M B and varying h A ■ M B would not necessarily 
give similar conductance variations. 43 

A. The quasiparticle current 

We now present numerical results on the quasiparticle 
current / qp given by equation (|30|) . In addition to the 
choice of Rashba spin-orbit coupling, we also need the 
angular dependence of the gaps X±,k- As before, we write 
the phase explicitly, such that x+,fc an d are real. 

We consider the same gaps as in Refs. @ andl20l given 
by r]k,s = ^ and r/k.T = A|n x fc|. The singlet and triplet 
components, rj^s an d %,Tj are defined in equations Q 
and . "J and A are treated as constants for simplicity. 
We also assume that "J > and A > without loss 
of generality. The gaps in the spin-orbit split bands are 
then 

X ±.k = *± A|n x fc|. (42) 

Let us define q = W/A. Whereas x+,k is fully gapped if 
q > 0, the gap X-,k contains line nodes if < q < 1. See 
Ref. [a for details. At this point, we should mention that 
other explanations of line nodes in CePtsSi have been 
put forwardi^ 

It is also for q < 1 that we may expect Andreev bound 
states at the surface, since x+,kX-,k < on a part of the 
Fermi surface in that case. However, one should note that 
formation of Andreev bound states does not depend on 
the presence of gap nodes. Isotropic x±.k with different 
signs will also result in subgap surface bound states. 

For simplicity, we assume a spherical Fermi surface and 
let the density of states be constant over the Fermi sur- 
face, N(k-p) = N. Let us introduce spherical coordinates 
by k = (cos (j>sm 8, sin sin 8, cos 8). As mentioned be- 
fore, the x-axis is perpendicular to the junction. In addi- 
tion, we let fi A and n B point along or opposite to the z- 
direction. The gaps are then given by X±,fe/A = (?±sin 8. 
For q < 1, Andreev bound states are formed for momenta 
with arcsin(g) < 8 < n — arcsin(g). As mentioned in 
section UlI A[ these surface bound states form below the 
smallest gap, i.e. below |x-.fc|/A = \q — sin#|. Figure[T] 
shows the spectrum of Andreev bound states E Q ^ in the 
case q = 0. The dependence on the asimuthal angle <fi 
is shown for three different polar angles 8. We see that 
E Q £ — > as 4> "~ * 5 which corresponds to k y = and 
thus |&fc|i | = 0. The maximal value of E Q f, is given by 
\q — sm9\. 

The tunneling matrix element, defined in equations 
and , will typically favour momenta with a large 
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FIG. 1: (Color online) Energy spectrum for the bound states 
E £ in the case q = 0. The dependence on the asimuthal 
angle (j> is shown for three different polar angles 6. 



component in the tunneling direction. Also, in the case 
of a smooth barrier, the parallel momentum is conserved 
in the tunneling process. Thus, we assume that 

\T k¥ , P A 2 = t kpxS(k\\ -P||) (43) 

will capture the qualitative features of the tunneling ma- 
trix element, where t is a constant 4^ 

One may show that the variation of the current with 
the angle £ disappears when only perpendicular momenta 
contribute. In other words, the effect is dependent on a fi- 
nite tunneling cone, where also nonzero parallel momenta 
contribute to the current. 

Tunneling spectroscopy on superconductors are inter- 
esting at low temperatures. At higher temperatures, 
sharp features giving information on pairing states may 
be smeared out. Thus, we investigate the limit of zero 
temperature here. However, for q < 1, the current at 
low voltage is dominated by resonant transport between 
Andreev bound states. This is contained in the sum 
g^pmv ' (kF,pF,iu} u ), where a product of two delta dis- 
tributions enters. At zero temperature, this leads to a 
discontinuity at V = 0, where the current jumps from 
zero to a finite value. The discontinuity disappears for 
nonzero temperatures and a sharp zero bias conductance 
peak appears. To get realistic current-voltage diagrams, 
we therefore retain a small temperature (T/A = 0.015) 
in this particular term, such that this discontinuity at 
zero voltage is smeared out. Such a small temperature 
will have no significant effect on the other terms. 

The current-voltage diagrams for several q are now pre- 
sented, where we have defined i qp = —I qp /(2eTT 2 t 2 N 2 ). 
We consider the cases of ( = and £ = it, i.e. equal and 
opposite directions of broken inversion symmetry. In ad- 
dition, we present the differential conductance G(eV) = 
d» qp /d(eV), which may be directly accessable in experi- 
ments. The latter has been obtained through a Savitzky- 
Golay smoothing filter 39 to remove noise from the numer- 
ical integration. 
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We start by considering the q = case, which corre- 
sponds to a pure spin triplet state. The gaps x+,fe an( i 
X-.fc are then of opposite signs on the entire Fermi sur- 
face except at fcp = ±n, where they have point nodes. 
Figure [5] shows the current-voltage diagram when q = 0. 
The differential conductance is presented in figure [H The 
large current at small voltages is due to transport be- 
tween Andreev bound states on each side. This gives 
rise to a zero bias conductance peak followed by nega- 
tive differential conductance. Similar phenomena appear 
in some <i-wave junctions.— We observe that there is no 
difference between the cases £ = and £ = ir in the pure 
triplet case. As stated before, this is also the case for the 
pure singlet case, q — > oo4^ However, we will see that 
this changes for finite q, when the gap is a mixture of 
singlet and triplet. 




eV/A 



2.5 



FIG. 2: (Color online) Current- voltage diagram in a Joseph- 
son junction when q — 9 / A — 0. Transport between Andreev 
bound states dominates for small voltages. There is no de- 
pendence on £ in this pure triplet case. 



70 

G(eV) 



35- 



? = 



0- 



q = 0.4. In this case, x+.k is fully gapped (although 
anisotropic) whereas the gap X-,fc has got line nodes at 
6 w 23.6° and 8 w 66.4°. Andreev bound states exist 
between these angles. Observe that the cases C = an d 
£ = 7r differs. This becomes clearer when studying the 
differential conductance in figure O We do not attempt 
to explain every feature in this figure, as this depends 
on the particular pairing state chosen. In addition, some 
of these features might also be smeared out in experi- 
mental results. However, the important thing to notice, 
which might be observable, is the qualitative difference 
of a junction with equal n-vectors (£ = 0) and one with 
opposite n-vectors (( = ir). 




eV/A 



FIG. 4: (Color online) Current-voltage diagram when q = 0.4. 
Transport between Andreev bound states dominates for small 
voltages. 
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FIG. 5: (Color online) Differential conductance in the case 
q = 0.4. Transport between Andreev bound states dominates 
for small voltages. 



FIG. 3: (Color online) Differential conductance as a function 
of voltage when q — 0. For small voltages, the zero bias 
conductance peak followed by negative differential resistance 
is due to transport between Andreev bound states on each 
side. There is no dependence on £ in this pure triplet case. 

Figure 3] shows the current- voltage diagram in the case 



The next current-voltage diagram, presented in figure 
|6l is for q — 1. Then, the line nodes have moved to the 
equator (6 — n/2) and will disappear for q > 1. Now, 
there is no part of the Fermi surface where x+,fcX-,fc < 0, 
such that there are no Andreev bound states. In the 
differential conductance in figure [71 there is a clear dif- 
ference between ( = and £ = it. See Ref. for a 
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simplified discussion of why this occurs. 



30r 

iqp/A 



22.5 




eV/A 



FIG. 6: (Color online) Current- voltage diagram when q = 1. 
At this point, there are no Andreev bound states. 
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G(eV) 



12.5 



■C = o 




eV/A 

FIG. 7: (Color online) Differential conductance when q = 1. 
Note the qualitative difference in the two cases for eV/ A < 4. 

Finally, we examine the scenario where the singlet to 
triplet ratio is q — 2. At this value, both x+,fc an d X—,k 
are fully gapped and of the same sign. The current- 
voltage diagram is given in figure [8] and the differential 
conductance in [5J Above eV/A — 2, the behaviour is 
similar to the q — 1 case. 

In the cases q — 1 and q — 2, we observe that the 
graphs differ in the region 2(q — 1) < eV/A < 2(q + 1). 
This will also be the case for higher values of q, but the 
width of this region (2(q + 1) — 2(q — 1) = 4) becomes 
small relative to the voltages at which the graphs differ 
(eV/A sa 2q). In the limit q — > oo, we are left with the 
singlet result, with a single step in the current for both 
( = and ( = n. 



B. The Josephson current 

We now move on to the Josephson current, given by 
equation (|37|) . This has not been investigated in as much 




FIG. 8: (Color online) Current-voltage diagram when q = 2. 
At this point, both bands are fully gapped. 
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FIG. 9: (Color online) Differential conductance when q — 2. 
The two cases of £ = and £ = n differ significantly when 
2 < eV/A < 6. 



detail as the quasiparticle current. In this section, we 
only suggest that the critical Josephson junction at zero 
voltage may depend on the angle ( between the axes of 
broken inversion symmetry of the crystal. We make no 
attempt to give any quantitative estimates here, since 
this depends not only on the particular pairing state of 
the material in question, but also on several other issues, 
like the details of the tunneling matrix elements. Only 
experiments can determine whether this effect really oc- 
curs and to what degree. 

The critical or maximal Josephson current at eV = 0, 
Ij,c(C) is defined as the absolute value of the Josephson 
current at phase difference $b — $A = ±tt/2. We still 
use the Rashba spin-orbit coupling and consider only one 
pairing state, given by 



X+m = const., x-,fe = 0. 



(44) 



This is probably not very realistic, but suffices to il- 
lustrate the effect 4& In this case, there are no Andreev 
bound states. 

As mentioned in the previous subsection, the depen- 
dence on C disappears when only perpendicular mo- 
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menta contribute to the current. This is also the case 
for the Josephson current. We illustrate this by intro- 
ducing a cutoff in the angle integrals, integrating over 
9 C < 9 < it - 9 C and -tt/2 + cf> c < (f> < n/2 - <fi c . Here, 
9 C = 4>c = corresponds to integration over the entire 
semisphere, whereas only perpendicular momenta con- 
tributes when 9 C — <f> c — > tt/2. 




°o —iii 



FIG. 10: (Color online) The variation of the critical Joseph- 
son current with C, for three different cutoff angles. Note the 
difference in the cases £ = and £ = 7r. The variation dimin- 
ishes as the tunneling cone is narrowed. 

Figure [TU] shows the variation of the critical Josephson 
current with Note the difference in current for the 
cases C = an d C = ^- One should also observe that 
the variation is reduced when the cutoff angle increases, 
corresponding to a narrowing of the tunneling cone. 

Although we have only studied a special scenario, the 
general message is that a variation of the critical Joseph- 
son current with £ may be expected when the gap is a 
mixture of singlet and triplet components. 

V. CONCLUDING REMARKS 

We have investigated both the current-voltage diagram 
and the critical Josephson current in planar tunnel junc- 
tions consisting of two superconductors with antisymmet- 
ric spin-orbit coupling. This is relevant for several re- 
cently discovered superconductors, where the spin-orbit 
coupling is a consequence of the crystal lacking inver- 
sion symmetry. Expressions for the currents have been 
derived in the tunneling limit using a general spin-orbit 
coupling. 

Numerical results have been presented in the case of 



the Rashba spin-orbit coupling. We have investigated 
the dependence on the relative angle between the direc- 
tions of broken inversion symmetry on each side of the 
junction. It has been shown that if the gap is a mixture 
of spin singlet and spin triplet parts, qualitative changes 
in the differential conductance may be expected when 
varying this angle. One may also observe quantitative 
changes in the critical Josephson current. This is a result 
of the fact that spin is conserved in the tunneling process 
whereas the spin structure of the spin-orbit split bands is 
determined by the direction of broken inversion symme- 
try. One should note that broken inversion symmetry on 
both sides of the junction is of importance. As stated ear- 
lier, similar conductance variations does not necessarily 
appear when replacing one of the superconductors with 
a ferromagnet and varying its magnetization. 

The experimental verification of these phenomena re- 
quire synthesis of junctions with specific crystallographic 
orientation on each side. It is worth mentioning that 
Josephson junctions with controllable crystallographic 
orientation were essential in proving the d-wave symme- 
try of the order parameter in the high-T c cuprates.— Fur- 
thermore, the roughness of the tunnel barrier should be 
as small as possible. In addition, the planar tunnel junc- 
tions must be thin enough to ensure that momenta with 
finite parallel components contribute to the current. Fi- 
nally, it should be pointed out that a difference in the 
normal phase densities of states of the two bands could 
give rise to some of the above mentioned effects even for 
conventional pairing. However, this should be possible to 
detect by measuring the current-voltage characteristics in 
the normal phase above T c . 

Many approximations and assumptions have been 
made in order to produce these results. Thus, the re- 
sults presented here are expected to be of qualitative 
value only. The main message is that experiments on 
Josephson junctions of noncentrosymmetric supercon- 
ductors may provide a direct connection between the 
possibly unconventional pairing and the lack of inversion 
symmetry in the crystal. 
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APPENDIX A: DERIVATION OF THE HALF-SPACE GREEN'S FUNCTION 

First, we note that the momentum space Gor'kov equations in the bulk are 



A(k,iw n )G h (k,iuj n ) = I, 



(Al) 
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where 

/ [ki)„ - (sk — //)] 1 — -Bfe • a -Afe \ 

A{k,iu n )=[ , (A2) 

V -AT []u n + (e k - fi)] 1 - B k ■ a* J 

and Gb(k,\u> n ) are matrices in spin times particle-hole space. The subscript b denotes bulk. uj„ = (2n + 1)tt/ f3 is a 
fermion Matsubara frequency. The definition of Gh(k,ioj n ) and the solution of equations (|A1[) are given in section 

nra 

We now want to determine the normal and anomalous Green's function when we restrict our system to a half- 
space, i.e x < 0. Contrary to the bulk case, the Green's function will not be diagonal in momentum space. We do 
however assume translational symmetry in the y- and z-directions, such that the Green's function will be diagonal in 
fey = k y y + k z z. It is convenient to work in a mixed basis, where we define the Green's functions Q aa i (x\ , x%, ku , r) = 
-( r TCsi,fcn,ff(7")4 a) . ,(0)) and J 7 aa ,(xi,x 2 ,k\\,T) = (T T c Xllfe||lCT (r) c X2> _ fe||lO -/(0)}. The Gor'kov equations in the 
continuum limit are 

o 

dxA(xi,x, —id x , k\\,vjj n ) G(x, X2,k\\,icj n ) = 5{x\ - x 2 ) 1 (A3) 

-oo 

in spin x particle-hole space. We have defined 

([iw„l - H N (x, -i9 x ,fe||)] 6(xi - x) -A(xx,x,k\\) 
i ( A4 ) 
-A^(x,xi,k\\) [hjj n l + H$,(x,-id x ,-k\\)] S{xi -x) 

where A and i?N are 2x2-matrices in spin space. The 4x4 Green's function is 

G(x, x 2 , fcj| , vjJ n ) -Fix, x 2 , fe|| , iu) n ) 
G(x,x 2 ,k\\,iw n ) = | _ | (A5) 

-T\x,x 2 ,k\\,\uj n ) -Q t (x 2 ,x, -fe||,-iw n ) 



and should fulfill proper boundary conditions. The equa- 
tions are valid for x±, x% < 0. The difference from the full 
space Gor'kov equations is the restriction x < in the 
integral. The bulk version of equation (|A3|) reduces to 

The pair potential A(xi , x, fen) in equation (|A3[) should 
be determined self-consistently. It is well known that it 
may differ from its bulk value near surfaces . 21 ' 22 However, 
we will now apply the usual approximation™ of replacing 
the pair potential by its bulk value. Although this is a 
crude approximation, it is expected to give qualitatively 
correct resultsi 17 ' 20 

One way of deriving the half-space Green's function 
is to consider an infinite system and then introduce a 
wall of infinitely strong nonmagnetic impurities in order 
to confine the electrons to one side of the system.— The 
wall of impurities must ensure that there is no trans- 
port ("hopping") across the wall and no interaction be- 



tween the two sides. Since we use a continuum model, 
a single plane of impurities at x — will provide an 
impenetrable surface. It will however not prevent inter- 
action between the two sides due to the possibly nonlo- 
cal nature of the pair potential. Nevertheless, this inter- 
action with "ghosts" on the other side of the impurity 
wall is tantamount to approximating A(xi,x,ku) by its 
bulk value. In other words, we construct an auxiliary 
system for x > such that a particle in x% "feels" the 
pair potential Ab(xi, x, ku) from all x as it would in the 
bulk. Thus, in the approximation stated above, we may 
extend the x-integral in equation (|A3|) to also include 
positive x and use the bulk pair potential Ab(xi, x, fen). 
However, we must demand that the boundary condition 
G(x\, X2, fell, = for X\ = or x 2 — is fulfilled, 
due to the infinitely strong impurities at x = 0. 

Having made the above mentioned approximation, it 
is easy to show that the ansatz 



G(xi,X2,k\\,fa) n ) = Gb(aJi,ar 2 ,fc||,ia; n ) - Gb(xi,0,fc||,ia; n )G b 1 (0,0, fe||,icj„)G b (0, a; 2 ,fc||,iw„) (A6) 

satisfies the boundary conditions and the Gor'kov equations. Thus, we have expressed the half-space Green's function 
in terms of bulk Green's functions. 
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Since we desire a description of the system in terms of plane wave states, we are interested in the Fourier represen- 
tation of the Green's function (IA5D. 



&(xi,X2,ku,iuj n ) 



dfc M / dAa^Gffci.fca.iwnJe- 1 * 1 --" 1 



(A7) 



Using the Fourier representation of the bulk Green's function, Gb(xi, ^2, i^ n ) = dk x &h(k, iw„)e lk *( Xl X2 \ 
we arrive at 



G(fe 1 ,fe 2 ,iwn) = G b (fci,iw„)(5(fci ) - G b (fei,iw„)(& b 1 (0,0, fc||,iw„)G b (A;2,iw„)| <S(fc 1; || - fc 2j ||, 



(A8) 



We see that the half-space Green's function differs from the bulk function by the second term, which is nondiagonal 
in the perpendicular components of the momenta. 



APPENDIX B: TRACE CALCULATIONS 



In section HlBl we defined cr. P 1 = a\ of a'i = 3. pl i t + a, p ^, ■ cr, where a\ = 1 + AJBfc • cr. Using the 

1*1,1*2 ■Ofc 1 fc || ^fc 2 K 1, K 2 1*1,1*2 t>k 

algebra of Pauli matrices, one arrives at 



fifX = 1 + MBt ■ < ) + M< -Bt 2 ) + -Bt 2 )+ iXpfiB* ■ (b* xB% 2 ), 



(Bl) 



and similarly for side B, where A — > B, A, p, /i — ► 7, 77, f and fc — > p. We now intend to find the functions Aj(fep,pp) 
defined in equation (|31[) . on which both the one-particle current I qp and the two-particle current I] depend. First, 
note that 



Tr rr A ~ — 9 

LT(T B£ a Pi,P2 - Z 

Tr a Xpp a< — 2 
p 

lr<7 fc 1 ,fc 2 (T pi,P2 Z 



1 + A 7 (££ • Bl 



"Pl ,P2 



f J k 1 ,k 2 l J Pi,P2 ' ^"fei.fca "pi,P2/ 



(B2) 



obtained by using Trl = 2 and Trui = 0. To simplify the expressions Aj(fep,pp), some useful relations are B^ f -b^ = 

6 F ' ^11 = ^fe I 2 anc ^ fcp Sf = ^l^fe I 2 — ^' ^ n Edition, we are only interested in Ai(k-p,p-p) as appearing in the 
Fermi surface integrals (|30|) and (|37[) . This allows further simplifications when using the symmetries XA.fc — X\k anc ^ 
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B 



-k 



Bk. Thus, when appearing in the integrals (|30|) and (|37p . we have 



A Xj (k F ,p F ) 
A x 2 ^(k F ,p F ) 
A^(k F , PF ) 

A xp ^(k F ,p F ) 



PF I 

1 + 7^ + »?(7 + v)\bl I 2 + A f 7 + v + 2 7 ^|^„ I 2 ) (BL • -B B 



^(1 - 71/)^^ • 6; 



Pl! 



1 + A/x + p(A + M )|b£ B I +7 (A + /i + 2Ap/x|6£ || | 2 (B£ F ■ ) + 7 p(l - XM ■ B 



PF 



1 + A/j + p(A + mM„ I 2 1 + iv + Vh + v)K t I 2 



1 - A/i + (2A/j + p(A + //)) |6£ y | 2 j [1 - 1V + {2 1V + 77(7 + 1/)) |b B | 2 
(A + p + 2A W i|tf ll | 2 )( 7 + ^ + 27^1^ | 2 ) + A 7 + / J (B£ f . ) 



+ [(A + M + 2A W |bfe ll 1^)77(1 ~ 7") + (A + M)r?(l + 7^)j (#£ F • 6* ) 
+ - A/x)( 7 + v + 27^|6« | 2 ) + p(l + Am) (7 + f)] (<, • Bp F ) 
+ [p(l - XfiMl - 71/) + p(l + A M )T7(1 + 71/)] (6^ • 6* ) 
+ (A^ + M7) (Bt-Bl) 



jp(X - n)n(j + yXfj, + ip(A + m) J [7^ + t>*7(7 + ") 



(B3) 



(BixBl).(BlxB 



B - 

PF 



The fermion Matsubara sums in 
the identity 



APPENDIX C: MATSUBARA SUMS 

and (138[) may be converted to contour integrals in the complex plane through 



1 



/3 



iw„)B(iw„) = -— <p dz A(z ~ \uj v )B{z)n F (z) 



(CI) 



for general A(z— \u v ) and S(^). The contour C must encircle the poles of the Fermi-Dirac function n F (z) = (l+e^ z ) _1 . 
The functions A(z — \oj v ) and B{z) appearing in section Um will have branch cuts and possibly poles on the lines 
Imz = \lo v and Imz = 0, respectively. This must be taken into account when deforming the contour. After the 
deformation has been performed, we may let iw„ — > eV + i0 + . 
The functions entering the sums (|32|) and (|38|) are 



| V B |2 

I A-y,p 



(z) = - 



'[ y B V B 4- V B V B — V B V B — 7 2 
'\ A7,pA^,p 1 Ar^pAz^p A'-y,pAl/,p ^ 



Ix^.J 2 - 



2(lx?,p| 2 - z) [b?, p|1 (x?, P X B , P - ^ 2 ) + 6?, p „ ^/lx?, P l 2 -^yixl^l 



■^■-vm/^) 



where we have defined c = sgn(7 + r] + v) and d = jrju. 
We choose the branch cuts such that 





B B B _ 2( B 1 V B _ B i 

A^y,pAr/,pAf ,p \A-y,p 1 Az^,p A?i,p 1 




V'lX-d.pl 2 - z 2 


2(lx?p| 2 -^ 




b B 

>pii 




+ ^P,Vlx?,pl 2 -^ 2 v / lx B ,pl 2 -^ 



'\ x * p \ 2 -(E±io+r 

and similarly for side A. 



|X B P | 2 - £ 2 6 |x B p| - |B| T isgn(£) ^ - | x B p |2 e \ E \ - \xl P 



(C2) 



(C3) 
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The functions g®(p,E + i0+) and T^{p, E + i0+) have the property B(E + i0+) = B(E - i0+)*. In addition, 



Imi?(_E + i0 + ) is an even function of E. Using this, one finds that the imaginary part of the sums in (f521) may be 
expressed as 



1m Si(eV + i0 4 



sgn(eU) 



AE Im 



A(E- \eV\ +i0+) Im B(E + i0+) n F (E - \eV\) - n F (E) . (C4) 

e 1 ^ , where 



The functions f*(p, E + iO + ) and A^p, E + i0+) have the property B(E ±i0+) = B K (E) ± iBj(E) 

the real functions Br(E) and B\{E) are even and odd in E, respectively. At eV = 0, this enables us to write the 
imaginary part of the sums in (I38p as 



ImSi(iO+) = 



sm(iT - i? A ) 



A\ poo 



dE 



A^B^E) + A I ( J B) J B R (^)j [l - 2n F (E) 
I 



(C5) 
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One might expect that this approximation breaks down 
when kp :X — + and that N± ih — > in that case, but this 
will not pose any problem here. 

There should also be an index A(B) on the cr's, since they 
depend on B^p)- We skip this, since the momenta (fc or 
p) tells us to which side they belong. 

In the case of a ferromagnet with magnetization M B , 
combinations like B£ F ■ M B + ■ M B = 2bt, ■ M B 
would not contribute when summing over all feii and us- 
ing b_k« = — &fei| ■ However, introducing some spatial 
anisotropy in the densities of states could change this. 
See also Ref. l40l where a junction between a 2DEG with 



Rashba spin-orbit coupling and a ferromagnet were con- 
sidered. 

The size of t could of course depend on £, at least if differ- 
ent junctions are used for different £. However, the quali- 
tive features of the current-voltage diagram would be un- 
affected by that. 

This statement depends on the two densities of states in the 
nonsuperconducting phase being equal. However, if they 
are not, this should give observable effects also above T c . 
These particular gaps have the nice property that the en- 
ergy integrals may be performed analytically, given some 
reasonable approximations. 



